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Abstract 

We study the first Kaluza-Klein excitation of the Higgs boson in universal extra dimensions as 
a dark matter candidate. The first-level Higgs boson could be the lightest Kaluza-Klein particle, 
which is stable due to the conservation of Kaluza-Klein parity, in non-minimal models where 
boundary localized terms modify the mass spectrum. We calculate the relic abundance and find 
that it agrees with the observed dark matter density if the mass of the first-level Higgs boson 
is slightly above 2 TeV, not considering coannihilations and assuming no relative mass splitting 
among the first-level Kaluza-Klein modes. In the case of coannihilations and a non-zero mass 
splitting, the mass of the first-level Higgs boson can range from 1 TeV to 4 TeV. We study also 
the prospects for detection of this dark matter candidate in direct as well as indirect detection 
experiments. Although the first-level Higgs boson is a typical weakly interacting massive particle, 
an observation in any of the conventional experiments is very challenging. 
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I. INTRODUCTION 



The particle identity of dark matter (DM) is one of the most important questions in 
particle physics, both from a theoretical as well as an experimental point of view. Since 
none of the particles in the Standard Model (SM) of particle physics can make up the 
DM, the problem points to new physics beyond the SM. The most popular class of DM 
candidates is weakly interacting massive particles (WIMPs), which are particles that have 
weak interactions and masses in the GeV to TeV range. 

One of the possibilities for new physics beyond the SM is provided by models with extra 
spatial dimensions. In particular, in models with universal extra dimensions (UED) 
the conservation of Kaluza-Klein (KK) parity ensures the stability of the lightest KK 
particle (LKP), which may be a viable WIMP DM candidate j^. In the five-dimensional 
UED model, the internal space has to be compactified on the orbifold S ,1 /Z 2 in order to 
give chiral fermions at the level of the zero modes and to avoid the existence of massless 
fifth components of the gauge fields. At the excited KK levels, each fermion field / in 
the SM is replaced by two Dirac fermions, /b and /s, having the quantum numbers of 
the corresponding left- and right-handed SM fermions, respectively. In addition, each of 
the KK excitations of the electroweak gauge bosons obtains a mass by eating a scalar 
that is a combination of the excitations of the SM Higgs field and the fifth components of 
the gauge bosons. In the limit of a large compactification scale R~ l in comparison to the 
Higgs vacuum expectation value, these Goldstone bosons are mainly composed of the fifth 
components of gauge bosons, and hence, the KK excitations of the charged component 
H ± and of the pseudo-scalar A are physical particles. 

The mass spectrum in the UED model, and thus the identity of the LKP, is affected 
by boundary localized terms (BLTs) in the Lagrangian, which are not determined by the 
model itself. In minimal UED (MUED) models, it is assumed that all such terms vanish 
at the cutoff scale A of the model, and are only generated at loop-level by renormalization 
group running. In the five-dimensional MUED model, the LKP is the first-level KK 
excitation of the U(l) gauge boson, the B 1 The phenomenology of this DM candidate 



has been extensively investigated in the literature [2l.l4l-ll2l]. See also the review in Ref. [13 

In models beyond the MUED model, where the assumption that the BLTs vanish 
at the cutoff scale is relaxed, the particle masses and interactions generally change. In 
Ref. 14j, it was shown, using a restricted set of BLTs, that the first-level Z boson, Z l , or 
the first-level neutral Higgs boson, H 1 , could be the LKP in such models. 

The phenomenology of Z 1 DM has been studied in the literature. The relic density and 
direct detection prospects were considered in Ref. jsj. Neutrinos from Z 1 annihilations in 
the Sun were studied in Refs. 15|, |l6[ and the gamma ray spectrum from Z 1 annihilations 



in Refs. 17], |l8j. In general, the detection of Z l DM is more challenging than for the B l , 



due to a larger preferred mass and a different distribution of annihilation channels. One 
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exception might be the gamma ray line signal, which receives large contributions from 
the gauge boson self-interactions. 

In this Letter, we study the remaining KKDM WIMP candidate, the H 1 . In addition, 
the pseudo-scalar A 0,1 has exactly the same properties as H 1 from the DM point of view. 
Hence, our results hold also for that DM candidate. The rest of the Letter is organized as 
follows: In Sec. [Til we calculate the relic abundance of H l DM, discussing also the effects 
of coannihilations. Then, in Sec. IHIt we estimate the direct detection prospects for this 
DM candidate. Next, in Sec. IIV| we consider indirect detection through the photon line 
signal and continuum spectrum, positrons, and neutrinos from DM annihilations in the 
Sun. Finally, in Sec. [Vj we summarize our results and state our conclusions. 



II. THE RELIC ABUNDANCE 

The standard calculation of the relic abundance of a thermally produced WIMP is 
reviewed, for example, in Ref. 19[. In the case that the mass splitting between the LKP 
and some of the other first-level KK modes is small, the effects of coannihilations are 



important 20[. Taking these effects into account, the abundance is given by 
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where x F = ttidm/Tf, moM is the mass of the DM particle, T F is the freeze-out tempera- 
ture, Mpi ~ 1.2 • 10 19 GeV is the Planck scale, and 
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is the effective number of relativistic degrees of freedom at freeze-out. The quantities 
and lb are given by 



I a = x F / a eS (x)x dx, 



h = 2x F I b c g(x)x dx. 
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Here, a e g and b e g are defined by the expansion cr e g = a cS + b cS v 2 + 0(v 4 ) of the effective 
cross section, 

= E^'^ 1 + A *) 3/2 ( 1 + A^V^+H (3) 

where o\j is the coannihilation cross section between the states i and j, gi is the number 
of degrees of freedom for the state i, and Aj = (m; — m LKP )/m LK p, with m LKP denoting 
the mass of the LKP. Finally, the freeze-out temperature is obtained from the relation 



Xp = In 



i 45 gm DM Mpi(a eS + 6b cS /x F ) 
2ir 3 ^g*x F 




(4) 



where g is the number of degrees of freedom for the DM particle, c ~ 1/2 is determined 
numerically, and 

fl* = J>(! + A*) 3726 "** (5) 



In this work, we have used the software package micrOMEGAs [21] to numerically calcu- 
late the relic density of H 1 DM, including coannihilations. In addition, we have checked 
the results analytically, using Eqs. ([I])-©, and found agreement between the two meth- 
ods. 

In the non-minimal UED model that we are studying, the mass spectrum has to be 
modified relative to the MUED model, so that H 1 is the LKP. Therefore, the rest of the 
first-level mass spectrum has to be fixed in some way in order to assess the effects of 
coannihilations, which depend strongly on the particle masses. Rather than considering 
the detailed effects of the BLTs, we take a more phenomenological approach by making 
an ansatz for the mass spectrum. Also, we do not consider modifications of the coupling 
constants due to the BLTs. The simplest ansatz for the mass spectrum is a universal mass 
splitting between the LKP and all the other first-level KK particles, parametrized by the 
relative mass splitting A = (mi — ^lkpV^lkp- This parametrization has previously 
been employed in the literature, e.g., in Refs. p, |7|. The resulting relic density, as well as 
the relic density calculated without coannihilations, is presented in Fig. [TJ We also show 
the 68 % confidence region Qh 2 = 0.1126 ± 0.0036 obtained from a combination of the 
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WMAP seven-year data, baryon acoustic oscillations, and supernovae observations 
If coannihilations are not important, the H l relic abundance falls within this region for 
mm — 2100 GeV. Coannihilations tend to increase this value, giving m#i — 2100 GeV 
for A = 10 %, mm ^ 2600 GeV for A = 3 %, and mm ^ 2700 GeV for A = 1 %. 

In Fig. [21 we show the separate effects of coannihilations with leptons and quarks, 
respectively. We find that coannihilations with leptons tend to decrease the preferred 
mass, while coannihilations with quarks tend to increase it. In the most extreme case 
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FIG. 1: The relic density of H 1 DM function of the mass m#i. The results are shown 
for the case without coannihilations as well as for coannihilations with a universal relative mass 
splitting A for all the first-level KK modes. 
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FIG. 2: The relic density of H 1 DM as a function of the mass m H i . Left panel: coannihilations 
with KK leptons only. Right panel: coannihilations with KK quarks only. The relative mass 
splitting for the coannihilating particles is given by A. 



that we consider, i.e., A = 1 %, lepton coannihilations might shift the mass to about 
1 TeV, while quark coannihilations could push it above 3 TeV. 
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FIG. 3: Tree-level Feynman diagrams contributing to H scattering on nuclei at quark level. 
The subscripts D and S on the KK quarks denote the SU(2) doublet and singlet Dirac fermions, 
respectively. 



III. DIRECT DETECTION 

Now, we turn to the experimental signatures of H 1 DM, starting with direct detec- 
tion experiments. The standard calculation procedure for WIMPs scattering on nuclei 
is reviewed in Ref. 
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In addition to the WIMP-quark scattering cross sections, the 
scattering rate depends on the WIMP distribution in the vicinity of Earth as well as on 
the structure of the nucleons in terms of quarks and gluons. Throughout this Letter, we 
assume a Navarro-Frenk-White (NFW) halo profile 24| with scale radius r§ = 20 kpc. 

The H 1 , being a scalar particle, scatters only spin-independently on nuclei. The tree- 
level Feynman diagrams contributing to H 1 -quaxk scattering are given in Fig. [3J The 
amplitude for the process H 1 (pi)q(p2) —> H 1 (ps)q(p4) is 



v 

M = -i^u 4 



(pi +P2)' 



2 
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m 1 



u 2 , 



(6) 



where y q is the Yukawa coupling of the quark flavor q to the Higgs field, Ui = u(pi), and 

m q i. In the non-relativist ic limit, p\ ~ ^3 ~ (m^i, 0) and 

1,2,3. Expanding the amplitude to lowest 



is 



we have assumed m n i 

u^ u 2 = 2m,^2, while u^Y u 2 = for i 
order in m q , we obtain 
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where 

Cq = Vq m q 

The WIMP-nucleus cross section is given by 
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where the reduced mass /i = m H im N / '{m H i +m N ), m N is the mass, Z the atomic number, 
and A the mass number of the nucleus, and 



E Cq pp,n | 2 Pi „ \ ^ C q 
777, i 77 ^— ' rn 



_q=u,d,: 



c,b,t q . 



(10) 



Here, = 0.020 ± 0.004, = 0.014 ± 0.003, f Td = 0.026 ± 0.005, = 0.036 ± 0.008, 



r 1( 



ft - ft ~ 0.84, and similarly /£ G ~ 0.83 [25 



= 0.118 ± 0.062, # G = 1 
Since the amplitude for the contribution from the quark flavor q is proportional to 
the square of the Yukawa coupling, oc (m q /v) 2 , where v is the vacuum expectation 
value of the Higgs field, the scattering is suppressed for all flavors except for the top 
quark. The heavy quarks q — c,b,t contribute to the scattering only through loop-level 
couplings to gluons. The effective couplings for these quarks, given in Eq. fllOp . were first 



derived in Ref. 26j] for neutralino-nucleon scattering, and they do not hold in general. 
Nevertheless, we use these expressions to estimate the contributions from the top quark, 
as the contributions from the light quarks are completely negligible for the H 1 . 

Experimental results are usually expressed in terms of the WIMP-nucleon cross section 
a n = a sl ml/(n 2 A 2 ). Using /| G ~ /£ G , we obtain 
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where A q = {m q i — m H i) /m^- For A q = 3 % and m H i = 2 TeV, the scattering cross 
section is several orders of magnitude below the sensitivities of current direct detection 
experiments, such as XENON100 27[, in the relevant mass range. The cross section 
increases with decreasing mass m^i and/or mass splitting A q . However, coannihilations 
with KK quarks drive the mass to larger values, as observed in Fig. [2j Hence, if A q is 
small, there is a tendency to drive the mass m#i to a value that is too large for successful 
direct detection. 
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Final state 


Branching ratio 


HH 


0.543 


ZZ 


0.237 


w+w- 


0.220 


ff 






TABLE I: Branching ratios into all final state channels for H DM. The branching ratios are 
computed in the limit of degenerate first-level KK masses. Annihilation into fermion-antifermion 
pairs is helicity suppressed, and thus, it is negligible due to the large compactification scale 
preferred for the relic abundance. 

IV. INDIRECT DETECTION 

Next, we turn to the indirect detection of H 1 DM through the observations of decay 
products from H 1 pair annihilations. We consider high-energy photons, positrons, and 
neutrinos from H l annihilations in the Sun. The model-dependent input to the indirect 
detection signals is the total annihilation cross section and the branching ratios into 
different final states. Since the typical WIMP velocity is v ~ 10~ 3 , these quantities are 
calculated at zero momentum, i.e., only s-wave contributions are taken into account. The 
H 1 annihilates dominantly into the final states HH, ZZ, and W + W~, and the total 
annihilation cross section is given by 

avi^H 1 ) ~ (0.83 pb) ( . (12) 

V m Hi ) 

For scalar DM, annihilations into fermion-antifermion pairs are helicity suppressed, i.e., 
the s-wave cross sections are proportional to (mf/mji) 2 . Since the relic abundance re- 
quires a relatively large value for R^ 1 , the annihilation cross sections are small even for 
the top quark. The branching ratios into all available SM final states are given in Table [H 
computed in the limit of a degenerate first-level KK mass spectrum. 

A. The gamma ray line signal 

A smoking-gun signature of DM would be the detection of a gamma ray line signal, 
coming from the loop-level process H l H l — > 77, with E 1 = mn\. In addition, the 
processes H X H X — > ^Z,^H could give line signals at E 1 = m H i[l — m 2 x j \Am 2 H1 )], where 
X = Z, H. The relative shift from the two-photon peak, AE 7 = — m 2 x / {Am 2 H1 ) ~ 10 -3 , 
is too small to be resolved experimentally, and hence, the individual peaks add up to 
a single one. For H l DM, s-wave annihilation into the ^H final state is not possible 
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FIG. 4: One-loop top quark mediated Feynman diagrams contributing to the process H H — > 
77. In addition to these diagrams, there are three diagrams with crossed final state particles. 
The subscripts D and S on the KK quarks denote the SU(2) doublet and singlet Dirac fermions, 
respectively. 



due to conservation of angular momentum, and therefore, the process is suppressed. The 
H l H l 7Z process is discussed below. 

The amplitude for the process H l {pi)H l {;p2) —> 7(^3)7(^4) can be written as 

M=M flu er( P3 )C(pi), (13) 

where e* are the photon polarization tensors. For annihilation at rest, p\ ~ P2 — p = 
(m#i, 0). Using conservation of 4-momentum, 2p + p 3 + p 4 = 0, the transversality of the 
polarization tensors, • e* = 0, and the Ward identity p^M^ = p iu M^ u = 0, the tensor 
J\A.^ V can be reduced to the simple form 

= b( V ^$- - IgV J , (14) 

where the quantity B depends on the particle masses only. The cross section is given by 

\B\ 2 

av = —. (15) 

In Fig. HI we show the Feynman diagrams for H l H l — > 77 that involve internal top 
quarks. From an analysis of the coupling constants only, we would expect these to be larger 
than the corresponding diagrams that involve internal bosons by a factor {yt/g) 4 — 5. 
From this subset of diagrams, we find that 

B = 8^(0, 0, 4, V , rj, rj) - 4C (1, 0, -1, 0, rj, rj) - 4, (16) 
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where rj = (m t i/m H i) 2 and 



6-0(0,0,4,77,77,77): 



arctan 



;i7) 



C (l, 0,-1, 0,77, 77): 
Here, Li 2 (x) denotes the dilogarithm, 

Li 2 (x) = 



— J - Li 2 ( - 



;is) 



1 log(l - xt) 



dt. 



(19) 



Finally, the flux at Earth in the direction of the galactic center is given by 28 



$ 7 ~ (4.7- 10~ 12 m" 2 s" 1 ) ( 



2 TeVV 

10~ 29 cm 3 s _1 ) V m H i J 



(Jgc)aqA^, (20) 



where the solid angle Afl represents the resolution of the detector and (Jgc)aq is the 
dimensionless line-of-sight integral in the direction of the galactic center. For AQ = 10 -5 , 
(Jgc) aqAQ ~ 0.13 for the NFW halo profile. In Fig. [51 we present the numerical results 
for a number of different values for the mass-splitting parameter 77. In addition, we 
have calculated the contribution from diagrams involving internal bosons numerically, 
and found a result which is of the same order of magnitude as that from the top quark 
diagrams. Hence, the order of magnitude of our result is correct. Due to the small size of 
the line signal, however, a more detailed calculation would not be useful. 

We have also studied the 7Z final state, and found that the corresponding photon flux 
is smaller than from the 77 final state by about one order of magnitude. This is partly 
due to the difference in couplings, and partly due to the fact that each H 1 !! 1 — > 7Z 
annihilation only results in a single photon. 



B. The photon continuum spectrum 



In addition to the line signal, there is a continuous component of the photon spectrum. 
Primary photons are produced in final state radiation (FSR) processes, i.e., three-body 
processes of the form H X H X — \ XX7, where X denotes an electrically charged SM particle. 
These processes are suppressed by a factor a ~ 1/137 relative to two-body annihilation 
processes. However, the addition of a photon in the final state opens up the possibility of 
annihilations into fermion-antifermion pairs, which are otherwise helicity suppressed. The 
situation is similar to the case of neutralino annihilations, where the Majorana nature of 
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FIG. 5: The integrated gamma ray flux from the process H 1 !! 1 — > 77 as a function of the LKP 
mass mfji, for different values of the mass splitting 77 = (m t i /m H i) 2 . 



the neutralino leads to helicity suppression. In addition to primary photons, secondary 
photons are produced in the decays of other final states, which are directly produced, 
e.g., quarks, leptons, and gauge bosons. Both of these contributions are calculated in 
micrOMEGAs. In Fig. [6j we present the total continuum spectrum as well as the separate 
contributions from primary and secondary photons. The hard part of the spectrum is 
dominated by the FSR. In contrast to the case of B l annihilations, there is no sharp 
cutoff at = m H i, due to the different distribution of annihilation products. The 
secondary photons generated by the decays of other particles mainly contribute to the 
soft end of the spectrum. 

The differential flux close to the endpoint E 1 = m H i is larger than the line signal by 
several orders of magnitude. Hence, the line signal is completely negligible in comparison 
to the continuum spectrum. 



C. Positrons 

In general, positrons from DM annihilations can be produced directly through the 
e + e~ annihilation channel as well as indirectly through the decays of other annihilation 
products. For the H 1 , the direct e + e~ channel is suppressed by the small Yukawa coupling 
constants for light fermions. Hence, the spectrum is dominated by secondary positrons, 
and therefore, it is relatively soft. 

The observed flux at Earth depends strongly on the propagation of the positrons 
through the galaxy. The most important effects are space diffusion and energy losses 



due to synchrotron radiation and inverse Compton scattering 29j. In micrOMEGAs, the 
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FIG. 6: The differential photon spectrum as a function of £L, for m H i = 2 TeV. 



propagation is modeled as a diffusion-loss equation for the number density of positrons 
per unit volume and energy, xjj e + = dn e +/dE, which is of the form 



Qe+(*,E). 



(21) 



Here, Q e + is the source term, the space diffusion coefficient K = Ko(E/Eo) ' 7 with Kq = 
0.0112 kpc 2 /Myr and E = 1 GeV, and the positron loss rate b(E) = E 2 /(E r E ) with the 
energy loss time te = 10 16 s. 

The positron flux at Earth is shown in Fig. [7] for m# i = 2 TeV. For comparison, we 
also give the corresponding results for annihilations of B l and Z 1 DM, for which we have 
chosen the typical masses m^i = 1 TeV and m^i = 2 TeV given by the relic abundance 
calculations for the respective LKP candidates. The B l spectrum has been shown to 
be possible to fit to the PAMELA data, although the predicted magnitude is too small, 



requiring a boost factor of about 10 3 30J. The H 1 spectrum is relatively soft and has no 
sharp cutoff at E e + = m H i. In comparison to the flux from B 1 annihilations, it is smaller 
by about two orders of magnitude. 



D. The neutrino spectrum 

The gravitational capture of WIMPs in the Sun could lead to a significantly enhanced 
density, giving neutrino signals from pair annihilations. However, the capture rate of 
WIMPs in the Sun is proportional to the WIMP-proton cross section, which is the same 
quantity that is constrained by direct detection experiments. This means that the ex- 
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FIG. 7: The differential positron spectra E 3 - 5 d$/(dE e +dn) for B 1 , Z 1 , and H 1 DM annihila- 
tions, with m^i = 2 TeV, m^i = 2 TeV, and m B i = 1 TeV, which are typical values for the 
respective LKPs if the correct DM abundance should be obtained. 



pected results in neutrino telescopes are correlated with the limits from direct detection 
experiments. An important consequence is that for DM candidates that interact only 
spin-independently with nuclei, the constraints from direct detection experiments are al- 
read y st rong enough to rule out an observation of neutrinos from DM annihilations in the 



Sun 31 . This is the case for scalar DM candidates, such as the H 1 , and hence, neutrinos 



from the Sun are not a promising detection channel for this DM candidate. 



V. SUMMARY AND CONCLUSIONS 

In this Letter, we have investigated the first-level KK excitation of the Higgs boson in 
non-minimal UED models as a DM candidate. We have calculated the relic abundance, 
including coannihilations with other first-level KK particles. In addition, we have studied 
the detection prospects in direct as well as indirect DM detection experiments. 

Although the H 1 is a typical WIMP DM candidate, we find that detecting it would 
be very challenging in direct as well as indirect detection experiments. The main reasons 
for the suppressed rates are the facts that the Yukawa couplings are small for all fermions 
except for the top quark and that annihilation into any fermion-antifermion pair is helicity 
suppressed. This means that the if 1 -quark coupling relevant for direct detection is small, 
that the continuum gamma ray spectrum does not display a sharp cutoff, and that the 
positron spectrum is soft, consisting mainly of secondary positrons. In addition, being a 
scalar, the H l interacts only spin-independently with nuclei, and is not expected to give 
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observable signatures in neutrinos from the Sun. 

To conclude, the DM phenomenology of the H 1 is very different from that of the 
standard KKDM candidate, the B 1 . The B 1 has relatively good detection prospects, 
especially in indirect detection experiments. The photon and positron channels both 
feature hard spectra with sharp cutoffs at -E 7 , e + = rUffi, and also, the mainly spin- 
dependent 5 1 -nucleon interactions give rise to potentially strong signatures in neutrinos 
from the Sun. The H 1 phenomenology is actually more similar to that of the Z l , which has 
a similar preferred mass range for the relic abundance, and also features large branching 
ratios into bosons, rather than fermions. The main exception is the gamma ray line signal, 
which is expected to be strong for the Z 1 , due to large contributions from the non-Abelian 
gauge boson self-interactions, but not for the H 1 . 

It is apparently extremely hard to positively identify the H 1 as a DM particle. However, 
it would in principle be possible to produce it at a collider like the LHC If collider 
experiments revealed part of the KK spectrum, with an if Mike particle as the LKP, while 
direct and indirect DM detection experiments only kept setting limits, then one would 
be very tempted to declare the H 1 as the DM particle. A more detailed experimental 
investigation of its properties could then be used to probe the agreement of abundance 
calculations with the observed value. In short, although the experimental situation is very 
challenging, there is still hope to be able to establish the H 1 as DM in the future. 

Throughout the Letter, we have assumed that the coupling constants are not affected 
by the BLTs. In general, this might not be the case, and the results could be changed by 
such effects. An investigation of this issue is beyond the scope of this work. 

Finally, we repeat that the results of this Letter can be directly carried over to the 
pseudo-scalar A 0,1 . 
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